CRITICAL AGE DEPENDENT BRANCHING MARKOV 
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Abstract. This paper studies: (i) the long time behaviour of the em- 
pirical distribution of age and normalised position of an age depen- 
dent critical branching Markov process conditioned on non- extinct ion; 
and (ii) the super-process limit of a sequence of age dependent critical 
branching Brownian motions. 

1. Introduction 

Consider an age dependent branching Markov process where i) each parti- 
cle lives for a random length of time and during its lifetime moves according 
to a Markov process and ii) upon its death it gives rise to a random number 
of offspring. We assume that the system is critical, i.e. the mean of the 
offspring distribution is one. 

We study three aspects of such a system. First, at time t, conditioned 
on non-extinction (as such systems die out w.p. 1) we consider a randomly 
chosen individual from the population. We show that asymptotically (as 
t — ► oo), the joint distribution of the position (appropriately scaled) and 
age (unsealed) of the randomly chosen individual decouples (See Theorem 
I2.ip . Second, it is shown that conditioned on non-extinction at time t, 
the empirical distribution of the age and the normalised position of the 
population converges as t — > oo in law to a random measure characterised 
by its moments (See Theorem 12. 2p . Thirdly, we establish a super-process 
limit of such branching Markov processes where the motion is Brownian (See 
Theorem E3]). 

The rest of the paper is organised as follows. In Section \2. II we define the 
branching Markov process precisely and in Section 12.21 we state the three 
main theorems of this paper and make some remarks on various possible 
generalisations of our results. 
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In Section [3] we prove four propositions on age-dependent Branching pro- 
cesses which are used in proving Theorem 12.11 (See Section [4]). In Section [3] 
we also show that the joint distribution of ancestoral times for a sample of 
k > 1 individuals chosen at random from the population at time t converges 
as t — > oo (See Theorem 13. 5p . This result is of independent interest and is a 
key tool that is needed in proving Theorem 12.21 (See Section [5]). 

In Section [6l we prove Theorem 12.41 the key idea being to scale the age 
and motion parameters differently. Given this, the proof uses standard tech- 
niques for such limits. Theorem 12.11 is used in establishing the limiting log- 
Laplace equation. Tightness of the underlying particle system is shown in 
Proposition 16.41 and the result follows by the method prescribed in [7]. 

2. Statement of Results 

2.1. The Model. 

Each particle in our system will have two parameters, age in R + and 
location in R. We begin with the description of the particle system. 

(i) Lifetime Distribution G(-): Let G(-) be a cumulative distribution 
function on [0, oo), with G(0) = 0. Let fi = sdG(s) < oo. 

(ii) Offspring Distribution p : Let p = {pk}k>o be a probability 
distribution such that po < 1, m = Yl'kLo ^Pk = 1 an d that a 2 = 
T,T=o k2 Pk - 1 < oo. 

(iii) Motion Process ??(•): Let rj(-) be a R valued Markov process 
starting at 0. 

Branching Markov Process (G,p,rj): Suppose we are given a reali- 
sation of an age-dependent branching process with offspring distribution p 
and lifetime distribution G (See Chapter IV of [5] for a detailed description). 
We construct a branching Markov process by allowing each individual to ex- 
ecute an independent copy of rj during its lifetime r starting from where its 
parent died. 

Let Nt be the number of particles alive at time t and 
(2.1) C t = {(4,Xj):i = l,2,...,N t } 

denote the age and position configuration of all the individuals alive at time 
t. Since m = 1 and G(0) = 0, there is no explosion in finite time (i.e. 
P(N t < oo) = 1) and consequently C% is well defined for each < t < oo 
(See E). 
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Let B(R + ) (and B(R)) be the Borel cr-algebra on R + (and R). Let M(R+ x 
R) be the space of finite Borel measures on R + x R equipped with the weak 
topology. Let M a (R+ x R) := {v G M(R+ xR):v = EtiK^i;-)^ G 
N, ai G R+, Xi G R}. For any set A G B(R+) and 5 G B(R), let y t (A x B) be 
the number of particles at time t whose age is in A and position is in B. As 
pointed out earlier, m < oo, G(0) = implies that Yt G M a (R+ x R) for all 
t > if Yq does so. Fix a function <fi G C^(R+ x R), (the set of all bounded, 
continuous and positive functions from R + x R to R+), and define 

, N t 

(2.2) {Y t A)= / ct>dY t = Y,m,Xl). 

J i=i 

Since rj(-) is a Markov process, it can be seen that {Yt : t > 0} is a Markov 
process and we shall call Y = {Y t : t > 0} the (G,p,rj)- branching Markov 
process. 

Note that Ct determines Y t and conversely. The Laplace functional of Y t , 
is given by 

(2.3) L t <P(a,x) := E a , x [ e -^ Yt >] = E[e~^ \ Y = 5 a , x }. 
From the independence intrinsic in {Yt : t > 0}, we have: 

(2.4) E ui+U2 [e~^] = (E Ul [e-^])(E U2 [e-^]), 

for any v { G M a (R + x R) where £7 1/ .[e - <*' y «>] := £;[e^^' ri ) | Y = n] for 
i = 1,2. This is usually referred to as the branching property of Y and can 
be used to define the process Y as the unique measure valued Markov process 
with state space M a (R + x R) satisfying L t + S (j)(a, x) = L t (L s (cj))){a, x) for all 
t,s > 0. 

2.2. The Results. 

In this section we describe the main results of the paper. Let A t be the 
event {N t > 0}, where N t is the number of particles alive at time t. As 
Po < 1, P(A t ) > for all < t < oo provided P(N = 0) / 1. 

Theorem 2.1. (Limiting behaviour of a randomly chosen particle) 

On the event A t = {N t > 0}, let (at,X t ) be the age and position of a 
randomly chosen particle from those alive at time t. Assume that rj(-) is 
such that for all < t < oo 

(2.5) E{n{t)) = 0,v(t) = E(r] 2 (t)) < oo, sup v(s) < oo, 

0<s<t 

poo 

and ip = v(s)G(ds) < oo. 
Jo 
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Then, conditioned on A t , (a*, ^|) converges as t —* oo, to ([/, V) m distri- 
bution, where U and V are Independent with U a strictly positive absolutely 
continuous random variable with density proportional to (1 — G(-)) and V is 
normally distributed with mean and variance ^. 

Next consider the scaled empirical measure G M a (M + x R) given by 
Y t (A xB) = Y t (A x VtB), A G B(R+), B G B(R). 

Theorem 2.2. (Empirical Measure) 

Assume A2.5\) . Then, conditioned on At = {Nt > 0}, the random measures 
{jf} converges as t — > oo in distribution to a random measure v, charac- 
terised by its moment sequence mk((f>) = E[v{cf)) k ], for cf> G C&(R + x M), 
k > 1. 

An explicit formula for m/%((/>) is given in (|5.2p below. 

Our third result is on the super-process limit. We consider a sequence of 
branching Markov processes (G n , p n , rj n ){ n>1 } denoted by {Y t n : t > 0}{ n>1 } 
satisfying the following: 

(a) Initial measure: For n > 1, Yg n = n nu, where ir nu is a Poisson 
random measure with intensity nv, for some v = ax fi G M(l + xM). 

(b) Lifetime G n (-): For all n > 1, G™ is an exponential distribution 
with mean j 

(c) Branching p n ,- ; F° r n > 1, Let F n (u) = YH?=oPn,kU k be the 
generating function of the offspring distribution p n = {p n ,k}k>o- 
We shall assume that F n satisfies, 

(2.6) lim sup || n (F n (l - u/n) - (1 - u/n)) - u 2 0, 
n ^°°0<u<N 

for all N > 0. 

(d) Motion Process rj n (-): For all n > 1, 

1 /"* 

(2.7) r? n (i = — / a(u)dB(u), t > 0, 

Vn Jo 

where {B(t) : t > 0} is a standard Brownian motion starting at 
and o : M + — ► R is a continuous function such that J °° o~ 2 (s)dG(s) < 
oo. It follows that for each n> 1, r] n satisfies (|2.5|) . 

Definition 2.3. Lei £ be an independent exponential random variable with 
mean ± ; < A < oo. For f G C, + (R + xl) letU t f(x) = E{f{£,x + ^fhpB t )) 
where ip is defined in \2. 5|) . For t > 0, Zei Ut(/) 6e the unique solution of 
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the non linear integral equation 

(2.8) u t f{x) = U t f(x) - A f U t - S {u s (ff ){x)ds. 

Jo 

Let {y*t ■ t > 0} be a M(R + x R) valued Markov process whose Laplace 
functional is given by 

(2.9) E Sxii [e-VM] = e-MM , 

where f G C, (R+ x R rf ) (^/ie sei of all continuous functions from R+ x R to 
R with finite limits as (a,x) — » ooj and Vj(/)(x) = Ut(f(x)) for i£l (See 
[7] /or existence ofy satisfying \2.9\) . 

Note that in the process : t > 0} defined above, the distribution of the 
age (i.e. the first coordinate) is deterministic. The spatial evolution behaves 
like that of a super-process where the motion of particles is like that of a 
Brownian motion with variance equal to the average variance of the age- 
dependent particle displacement over its lifetime. Also, u s (f) in second 
term of (|2.8p is interpreted in the natural way as a function on R + x R with 
u s (f)(a, x) = u s (f)(x) for all a > 0, x G R. 

Theorem 2.4. (Age Structured Super-process) 

Let e > 0. Let {Y t n : t > 0} 6e the sequence of branching Markov processes 
defined abovefi.e.in (a), (b), (c), (d)). Then as n — > oo ; {y™ = ^Y™ t ,t > e} 
converges weakly on the Skorokhod space D([e, oo), M(R + x R)) to {3^ : t > 

2.3. Remarks. 

(a) If r/(-) is not Markov then C t = {a l t , Xl,fjt,i = {^t,i(^) : < u < 
a\} : i = 1,2... , iVf} is a Markov process where {fjt,i(u) ■ < u < a\} is 
the history of rj(-) of the individual i during its lifetime. Theorem 12.11 and 
Theorem 12.21 extends to this case. 

(b) Most of the above results also carry over to the case when the motion 
process is R d valued (d > 1) or is Polish space valued and where the offspring 
distribution is age-dependent. 

(c) Theorem 12.11 and Theorer d2.2l can also be extended to the case when 
r](Li), with L\ = G, is in the domain of attraction of a stable law of index 
< a < 2. 
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(d) In Theorem 12.41 the convergence should hold on D([0, oo),M (M+ x R)) 
if we take a in the sequence of branching Markov processes to be £ (i.e. 
Exponential with mean j). 

(e) The super-process limit obtained in Theorem 12.41 has been considered 
in two special cases in the literature. One is in [6] where an age-dependent 
Branching process is rescaled (i.e. the particles do not perform any motion). 
The other is in [8] where a general non-local super-process limit is obtained 
when the offspring distribution is given by p\ = 1. In our results, to obtain 
a super-process limit the age-parameter is scaled differently when compared 
to the motion parameter giving us an age-structured super-process. 

(f) Limit theorems for critical branching Markov processes where the mo- 
tion depends on the age does not seem to have been considered in the liter- 
ature before. 

3. Results on Branching Processes 

Let {Nt : t > 0} be an age-dependent branching process with offspring 
distribution {pk}k>o and lifetime distribution G (see [5] for detailed discus- 
sion). Let {Cfc}fc>o be the embedded discrete time Galton- Watson branching 
process with fa being the size of the fcth generation, k > 0. Let At be the 
event {N t > 0}. On this event, choose an individual uniformly from those 
alive at time t. Let M t be the generation number and a t be the age of this 
individual. 

Proposition 3.1. Let At,at,Mt and Nt be as above. Let fi and a be as in 
Section \2.1\ Then 



2 



(a) lim tP(A t ) = 

t— >oo (J 



Nt 2 t ix 

(b) For all x > 0, lim P{— > x\A t ) = e ~Z r , 

t— »oo t 

M 1 

(c) Foralle>0, lim - -| > e\A t ) = 

t^oo t fl 

(d) For all x > 0, lim P(a t < x\A t ) = - f (1 - G(s))ds. 
Proof : For (a) and (b) see chapter 4 in [5]. For (c) see [9] and for (d) see 

m- □ 

Proposition 3.2. (Law of large numbers) Let e > be given. For the 
randomly chosen individual at time t, let {Lu : 1 < i < Mt}, be the lifetimes 
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of its ancestors. Let h : [0, oo) — > M be Borel measurable and E(\ h(L\) |) < 
oo with L\ = G. Then, as t —>■ oo 

-. M t 

p{ ~\ w E h ^ - E ^))\ > e ^*) - °- 

' 8=1 

Proof : Let e and e\ > be given and let fei(i) = — e) and k 2 {t) = 
£(i + e). By Proposition 13.11 there exists 5 > 0, rj > and to > such that 
for all t > to, 

(3.1) iP(iVi > 0) > 5 and P(N t < tn\A t ) < e x ; 



(3.2) P{M t e[k 1 {t),k 2 {t)f\A t )<e l . 

Also by the law of large numbers for any {Lj}j>i i.i.d. G with E\h(L\)\ < oo 

1 j 

(3.3) lim P(sup -| V h{Li) - E(h{L x ))\ > e) = 0. 

Let {Ck}k>o be the embedded Galton- Watson process. For each t > and 
A; > 1 let 0^ denote the number of lines of descent in the k-th generation 
alive at time t (i.e. the successive life times {-Li}i>i of the individuals in 
that line of descent satisfying J2i=i Li < t < J2i=i Li)- Denote the lines of 
descent of these individuals by {Cktj '■ 1 < j < Cfct}- Call Cktj bad if 

k 

(3.4) \~J2h(L ktji )-E(h(L 1 )))\>e, 

i=i 

where {Lktji}i>\ are the successive lifetimes in the line of descent Cktj start- 
ing from the ancestor. Let Ckt,b denote the cardinality of the set {Cktj '■ 1 < 
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j < Ckt and Cktj is bad}. Now, 



, M t 

p(\ w Y. h ^)- E ^))\>^) 



'M t f , 
i=i 



P( The chosen line of descent at time t is bad \A t ) 



< P( The chosen line of descent at time t is bad, M t £ \k\{t), k2(t)])\A t ) 
+P(M t e [ki(t)),k 2 (t)] c \A t ) 

1 2^j=k 1 (t)^jt,b 



E{ J ~Z' ;A t ) + P(M t e Mt)), k 2 (t)] c \A t ) 



P{N t >0) v N t 
" p(iv t >o)^ ( ivj ^ t >t V ) + 

<yk 2 (t) > 

+ m>0) ^( J= ^ } ;^ < */) + P(M t e [*!(*)), fe(t)] e |A*) 

+ P p{N t ~>0) + P{Mt ^ ** 
1 fa(*) 

= ^ > o) . + 

+P(iV t < ty|# t > 0)+P(M t e [A;i(t)),fc 2 (i)] c |^ t ) 

(3.5) 



For t > to by f|3.2|) and ()3.3j) . the last two terms in (|3.5p are less than e\. 
The first term is equal to 



1 fca(t) . fca(t) 

ir?P(iV t > 0) S = trjP ( Nt > o) . S ^E 1 ^ is bad.}) 

j—k±(t) j—k\(t) i—1 
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1 fca(t) 

tr,P(N t > 0) ,.^/ (Ci) 



xP ^ Li < * < L *> "I S ^ " ^MLi))! > e J , 
\i=l i=l ^ i=l / 

where the {Lj}j>i are i.i.d. G. 
Using (|3.ip and (since m = 1) E(Qj) = i?(Co) we can conclude that 

, te(t) 

tr)P{Nt > o) £ 

/ ^ v p ( BU Pi>M*)jlS*=iM^)-^(M^i))l >e) 

< -MCo) 



*7/P(iV t > 0) 

P(SUP„^ m i| 

< £(Co) 



P ( su P.>fciW jl H=i ^(Li) - £(MLi))| > e) 



7/5 

(3.6) 

which by (|3.3|) goes to zero. So we have shown that for t > to, 



1 Mt 



i=l 

Since t\ > is arbitrary, the proof is complete. □ 

Proposition 3.3. Assume $2.5]) holds. Let {Lj}j>i be i.i.d G and {r/i}i>i 
be i.i.d copies of r\ and independent of the {Lj}j>i. For 6 £ M, t > define 
4>(0,t) = Ee 1 ®^. Then there exists an event D, with P(D) = 1 and on D 
for all 0ER, 

j^(pl—=,Ljj —> e 2 , as n — > oo, 



. jn 



where ip is as in \2. 5\) . 

Proof: Recall from (j23j) that u(i) = E(r/ 2 (t)) for t > 0. Consider 



X ni = . Vi ^ Li/ lor ! < / < „ 



^l^i) 

and ^ = o"(Lj : i > 1). Given J 7 , {X n { : 1 < i < n} is a triangular array 
of independent random variables such that for 1 < i < n, E{X n i\T) = 0, 
T.UE{Xl l \T) = l. 
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Let e > be given. Let 



L n {e) = Y J E{X 2 ni ;X 2 ni >e\F) 



i=l 

By the strong law of large numbers, 

(3.7) — > ip w.p. 1. 

n 

Let D be the event on which (13.71) holds. Then on D 



limsup L n (e) < limsup V E^iU)^) :| Vl {Li)\ 2 > 



< limsup %E(]rn (Li)\ 2 :| 771C&1) | 2 > fc) 
= 0. 

Thus the Linderberg- Feller Central Limit Theorem (see [I]) implies, that on 
D, for all 9 € R 



i=i Vv^ivC^i) 



2 



Combining this with f|3.T|) yields the result. □ 

Proposition 3.4. For i/ie randomly chosen individual at time t, let 
{Lu,{rju(u) : < u < Lu} : 1 < i < Mt}, be the lifetimes and motion 
processes of its ancestors. Let Z t \ = -7= Yli^i Vti(Lti), and 
C t = a{M u L ti :l<i< M t }. Then 

(3.8) E (\E{e ieZtl \C t ) - e"T* | \A^j -► 

Proof: Fix 6 E K,ei > and e > 0. Replace the definition of "bad" in 
433} by 



(3-9) I Yl <K-t=, Lktji) ~ e"^| > e 

i=i vfc 



By Proposition 13.31 we have, 



(3.10) lim P(sup|TT0(^=,Li) -e"~| > e) = 0. 
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Using this in place of (13.30 and imitating the proof of Proposition 13.21 
(since the details mirror that proof we avoid repeating them here), we obtain 
that for t sufficiently large 

(3.11) P(| J]0( L U ) - e-T*! > eMt) < e. 

i=l * * 

Now for all flel, 

E(e^\£ t ) = l\0(^=,Lu). 

i=i * * 

So, 

limsup£(|£(e v^t 2 -^ tlj |£ t ) - e 2 ||A t ) 

= limsup£(| TT </>(^=, L«) -e-V||A t ) 

Aft a g2 
< ei + 21imsupP(|TT(/>(^==,L 4l )-e~^ i | > <n\A t ) 

2=1 

= ei + 2e. 

Since e > 0, ei > are arbitrary we have the result. □ 
The above four Propositions will be used in the proof of Theorem 12.11 

For the proof of Theorem 12.21 we will need a result on coalescing times of 

the lines of descent. 

Fix k > 2. On the event A t = {N t > 0}, pick k individuals C%, C2, . . . , Ck 

from those alive at time t by simple random sampling without replacement. 

For any two particles Cj, Cj, let rc^c^t be the birth time of their most 

recent common ancestor. Let Tk-x,t = sup{rcg,Cj,t '■ i ^ j, 1 < i,j < k}. 

Thus Tk—it is the first time there are k — 1 ancestors of the k individuals 

C\ , C2, ■ ■ ■ , Cf. . More generally, for 1 < j < k — 1 let r 3 - i as the first time 

there are j ancestors of the k individuals C\ , Ci , ■ ■ ■ Cfc. 

Theorem 3.5. 

(i) For any lim^oo P\ — 1 -j 21 - < x\A t ) = H(x) exists for all x > 
and H(-) is an absolutely continuous distribution function on [0, 00] 

(ii) Conditioned on At the vector ft = j(Tjt :l<j<fc — 1) as t — > 00 
converges in distribution to a random vector T = (Ti, . . . , Ifc_x) 
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with < T\ < T2 < . . . < Tjfc_i < 1 and having an absolutely 
continuous distribution on [0, 1] ^ 

Proof : The proof of (i) and (ii) for cases k = 2, 3 is in [9]. The following 
is an outline of a proof of (ii) for the case k > 3 (for a detailed proof see [3]). 

Below, for 1 < j < k — 1, Tj t t will be denoted by Tj. It can be shown that 
it suffices to show that for any 1 < i\ < i<i . . . < i p < k and < r\ < r^ < 
. . . < r p < 1, 

lim P(^ < ri <^<r 2 <...<^<r p <^i< r k . x < l\A t ) 

t-*co t t t t 

exists. We shall now condition on the population size at time tr\. Suppose 
that at time tr\ there are n\\ particles of which k\\ have descendants that 
survive till time tr%. For each 1 < j < k\\, suppose there are n 2 j descendants 
alive at time tr 2 and for each such j, let k 2 j out of the n 2 j have descendants 
that survive till time tr%. Let k 2 = (&21) ■ ■ ■ > &2|fci|) an d \k 2 \ = Ylj=i^2j- 
Inductively at time trj, there descendants for the j-th particle, 1 < 

j < |&£-i|. For each such j, let k{j out of riij have descendants that survive 
up till time tri + \ (See Figure [3] for an illustration). 
It will be useful to use the following notation: Let 

Tin, &n G N, fcn < nn, \ki \= ku,m = (nu). 

For i = 2, . . . i p let (m, h) G N;, where N t = N^- 1 ' x N^ 1 " 1 ' 

hj < n-ij, I h 1= hjy ( /- ) — II ( fc 3 . ) ■ 

j=l \ */ j=i \ w 



Let / a = P(iV s > 0). Now, 



P(f<n<f<r 2 <...<^<r p <^< r fc _i < 



^(n^^W ; /in 

II ff (?) (/^) fcy (! - UT^mL = mjWL > 0) x 



i=l 3=1 

rfc 



x^(k)^^|— , 

with ttj = r i+ i — rj, i = 1, 2, ... ,p — 1, u p = 1 — r p , iVj? u . is number of 
particles alive at time tu{ of the age-dependent branching process starting 
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with one particle namely j, g(k) = g(k\, . . . , k p ) is the proportion of con- 
figurations that have the desired number of ancestors corresponding to the 
given event, X? I N^N^ > and S = E^ 1 ' xj ■ 

Let qi = for 1 < i < p — 1. Then following [9] and using Proposition 
13.11 fi). (ii) repeatedly we can show that P(^- < T\ < ^f- < T2 < ■ ■ ■ < -j- < 
r p < ^j^- < Tk-\ < t\A t ) converges to 



e~ xqi x 



p+1 |fei-l| . , , k .j 

II II / dxe~ A -^e-^)g(k) 



x 

i=2 j=l 



8=1 



1 p+1 (<?j)' fcil 

- e n (ir^wn^d 5(k) x 

, fc+1 

n ^ x 



i=l 



Consequently, we have shown that the random vector ft converges in 
distribution to a random vector T. From the above limiting quantity, one 
can show that the T has an absolutely continuous distribution on [0, 1] . 
See [3] for a detailed proof. 

□ 



4. Proof of Theorem 12.11 

For the individual chosen, let (at,X t ) be the age and position at time t. 
As in Proposition 13.41 let {L t i, {r]u(u), < u < L t i} : 1 < i < M t }, be the 
lifetimes and the motion processes of the ancestors of this individual and 
{ r 1t(M t +i)( u ) : < u < t — J2i=i^ti} be the motion this individual. Let 
Ct = a(Mt, Lu,l < i < Mt). It is immediate from the construction of the 
process that: 

Mt 

at = t — L t i , 

i=l 
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Figure 1. Tracking particles surviving at various times 



whenever M t > and is equal to a + 1 otherwise; and that 

M t 

X t = X + ^2 Vti{L t i) + Vt(M t +i)(at)- 
i=i 

Rearranging the terms, we obtain 

«*" 7i ] = (a " \fl z " } + (0 ' (vr - V^) Zn) + <0 ' Tt + Za) ' 

where Z n = ^ — and Z t2 = -^Vt(M t +i)( a t)- Let e > be given. 
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P{\Z t2 \ > e\A t ) < P{\Z t2 \>e,a t <k\A t ) + P{\Z t2 \>e,a t >k\A t ) 
< P(\Z t2 \> e,at<k\A t ) + P(a t > k\A t ) 



< mZt2 ^ { a ^ k]At) +P(a t >k\A 



e 



By Proposition 13.11 and the ensuing tightness, for any r/ > there is a k^ 

P(a t > k\A t ) < |. 

for all k > k v ,t > 0. Next, 



E(\Z t2 \ 2 I at < kv \A t ) = E(I at < kv E(\Z t2 \ 2 \C t )\A t ) 
= E(I at < kri —^-\A t ) 
sup u <fc„ v(u) 

t 

Hence, 

P(\Z t2 \ > e\A t ) < lu fj +\ 
Since e > and rj > are arbitrary this shows that as t — > oo 
(4.1) Z t2 |A t 0, 

Now, for X > 0,9 £M, as at is £j measurable we have 

E{e~ Xat e~ l ^ Ztl \A t ) = E(e~ Xat (E(eT iez ^ \C t ) - e~^)\A t ) + 

e 2 i> . 

+ - — E( e - Xat \A t ) 



Proposition 13.31 shows that the first term above converges to zero and 
using Proposition 13.11 we can conclude that as t — > oo 

(4.2) (at ,4= Z tl )\A t (U,V) 

v ^ 

As Xo is a constant, by Proposition 13.11 (c). (|4.2[) . (j4.1j) and Slutsky's The- 
orem, the proof is complete. □ 
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5. Proof of Theorem 12.21 

Let (i 6 Cj(K x We shall show, for each k > 1, that the moment- 

functions of g( <y \A t ) converges as t — > oo. Then by Theorem 16.16 in 
[TT] the result follows. 

The case k = 1 follows from Theorem 12.11 and the bounded convergence 
theorem. We shall next consider the case k = 2. Pick two individuals C\, C% 
at random (i.e. by simple random sampling without replacement) from those 
alive at time t. Let the age and position of the two individuals be denoted 
by (a\,Xl),i = 1,2. Let r t = Tc lt c 2 ,t be the birth time of their common 
ancestor, say D, whose position we denote by X Tt . Let the net displacement 
of C\ and C2 from D be denoted by X\_ TV % = 1,2 respectively. Then 
Xi = X Tt +Xl n ,i = l,2. 

Next, conditioned on this history up to the birth of D{= Qt), the ran- 
dom variables (a\,Xl_ n ),i = 1,2 are independent. By Proposition 13.51 (i) 
conditioned on At, converges in distribution to an absolutely continuous 
random variable T (say) in [0,1]. Also by Theorem 12.11 conditioned on Q% 

X* 

and A t , {{a\, /~ Tt ) > ^ = lj 2} converges in distribution to {{Ui, Vi),i = 1, 2} 

which are i.i.d. with distribution (U, V) as in Theorem 12. 11 Also -7= condi- 
tioned on A Tt converges in distribution to a random variable S distributed 
as V. 

x^ 

Combining these one can conclude that {{a\, -y=), i = 1,2} conditioned on 
A t converges in distribution to {(Ui,VTS + — T)Vi),i = 1,2} where 
U\,U2, S, V\, V2 are all independent. Thus for any <f> G C&(R + x M) we have, 
by the bounded convergence theorem, 

(5.1) 




Using Proposition 13.11 (b) and the fact that (f> is bounded we have 
lim^oo E({ Yt ^ ) 2 \A t ) exists in (0, 00) and equals 777-2(0). The case k > 2 
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can be proved in a similar manner but we use Theorem 13.51 (ii) as outlined 
below. First we observe that as <fi is bounded, 



E 



<Yt N i >k \A t J+ = Y^Eh(N t ,k) ( 11^,^)1 A t I +g(^C t ,N t ), 



where h(N t , k) — > 1 and g((p, C t , N t ) — > as t — > oo; and i = {ii, • • • > ^fc} 
is the index of ^particles sampled without replacement from Ct (see (|2.ip ). 
Consider one such sample, and re-trace the genealogical tree T\ £ T(k),(T(k) 
is the collection of all possible trees with k leaves given by i), until their most 
common ancestor. For any leaf ij in Ti, let 1 = n(ij, 1) < n(ij, 2) < • • • < 
n^jjNij) be the labels of the internal nodes on the path from leaf ij to 
the root. We list the ancestoral times on this by {tl, r n ^ i^, . . . ,T n ^^ Ni y 
Finally we denote the net displacement of the ancestors in the time intervals 

[0, Ti], [Ti,T n (i jt 2)], ■■■ , [T"n(ij,Ar i .-l)) T n(ij ) Ar i .)]' i T n(ij M ■ ) > *] 

by 

Vi, (n ) , fji. (r n(ij >2 ) , ri ) , . . . , ^ Ij (r n(ij ^ ) , r n(ij ^ _i) ) , fj-. (t, r n{ij )Ar . . } ) . 
Given the above notation we have: 



where 



1 

f(0J,t) = <j){a t 3 .—{fjl^Ti) + E ^(7H(< J -,m)»^(< J -,m-l)) + ^ (*> r n(*j.^))- 



m=2 



Now by Theorem 13.51 

) T n(i,.N; .)) d 

At — ► {Ti , T n ^. 2 ) , ■ ■ ■ , T n < yi] > AT i . ) ) • 



( r l > T n(ij,2) ) ■ • • j T n(ij ,N ij ) ) 



Vt 

So by Theorem 12.11 



lim^j^M^j |A t ) = J B(EEII^»i>*)l^l 



i Te7li=i 



(5.2) 
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where 

U,S\/Ti+ ^ ^J^(ij,m) - ?n(ij,m-l) + a/"*- — ^n(ij,Ni) 



l 3 



m=2 

with S, Z^., Z™, m = 2,... ,JV^, are i.i.d.l/, {/ is an independent random 
variable given in Theorem 12.11 and Tj's are as in Theorem 13.51 (ii). Since 
4> is bounded, the sequence {rrik(4>) = hmt->oo E ( <Y Vt > ) } is necessarily a 
moment sequence of a probability distribution on K. This being true for 
each <p, by Theorem 16.16 in [TT] we are done. □ 



6. Proof of Theorem 12.41 

Let Z be the Branching Markov process Y described earlier, with lifetime 
G exponential with mean A, p\ = 1 and r\ = r\\ (see (|2.7p ). Then it is easy to 
see that for any bounded continuous function, St4>(a, x) = E/ a x ) < Z t , 4> >= 
E(a,x)4'( a t, X t ) satisfies the following equation: 

-><txTr ±r „ _\ , / j„\ — As 



(6.1) S t <j>(a,x) =e- M W t <f>(a,x)+ [ ds\ e - As W s (S^ s (4>){0,-))(a,x), 

Jo 

where Wt is the semi-group associated to rjx. Let C be the generator of 
7/i . Making a change of variable s — > t — s in the second term of the above 
and then differentiating it with respect to t, we have 

j t S t {(t>){a,x) = -\e- xt W t 4>{a,x) + e~ xt CW t 4>{a,x) + XS t ((f>)(0,x) 



[ daX(-Xe- x ^)W t -.(S a (<f>)(0,'))(a,x) 
J o 



'0 

+ f d S Ae- A (*- s )£m_ s (,S s (0)(O,-))(a,x) 
Jo 



XS t (<j>)(0,x) 
+ (C-X) 



er xt W t cP(a,x)+ / ds\e~ x ^W t „ s (S s (cP)(0,-))(a,x) 



\S t (<t>)(0,x) + (£ - \)S t (<f>)(a,x) 



^l(a,x) + ^MAS^(a,x)+X(SM(0,x) - S t (<f>)(a,x)), 
oa 2 



AGE DEPENDENT BRANCHING MARKOV PROCESSES 19 

For each n > 1 define (another semigroup) K^4>(a, x) = E at o((j)(at,x+^%). 
Now note that, 

X 

Rt4>(a,x) = E a0 (4>(a t , x + —= ) 

V n 



in 

= S t <pn(a, y/nx), 

where 4> n {a,x) = <p(a, -j?=)- Differentiating w.r.t. t, we have that the gener- 
ator of i?" is 

(6.2) K n ^{a, x) = ^(a, x) + ^^A0(a, x) + A(0(O, x) - 0(o, x)). 

Proposition 6.1. Let e > and £ > e. Let (p e C Z + (R + x R d ). T/ien, 

(6.3) sup | i?^(0)(a,x) - U t (4>)(x) |-> 0. 

(a,i)el + xl 

Proof: Let t > e. Applying Theorem 12.11 to the process we have 
(a n f, ^=0 — (C^, V). The proposition is then immediate from the bounded 
convergence theorem and the fact that <fi G C ; + (M+ xl) □ 

Proposition 6.2. Let ir nu be a Poisson random measure with intensity nv 
and t > 0. The log-Laplace functional ofy™, 

(6.4) E nnv [e-(*W) = e-W\ 
where 

(6.5) u^{a,x) = R n nt n(l-e-^){a,x)-\ f dsR n n{t _ s) {nH n C^)){a, x), 



o 



where 

¥„(0)(a, x) := [F n (l - 0(0, x)) - (1 - 0(0, x))] . 

Proof: For any n € N, let Y^ 1 be the sequence of branching Markov pro- 
cesses defined in Section \2.2l It can be shown that its log-Laplace functional 
Lf satisfies, 
(6.6) 

rnt 

L^(a,x) = e- Xnt W^ t [e-' t> ](a,x)+ / dsXe~ Xs W^ [F n (L n nt _J(0, •))] (a,x)ds, 

Jo 

where t > and W^ 1 is the semigroup associated with rj n . Using the fact 
that e~ Xu = 1 — Jq U ds\e~ Xs for all ti > and a routine simplification, as 
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done in [10J, will imply that 
(6.7) 

nnt 

L^(a,x)=W^ t [e-^](a,x) + X / tC_ s (F n (L>(0, •)) - L»(a, x)ds 



o 



Therefore v™ t (<j))(a, x) = 1 — L™cj)(a,x), satisfies, 
(6.8) 

pnt 

v^(a,x) = WS(l-e-*)(a,x)+ / d S W%- a ((l-^0)-^«(l-«?)0(O, 

JO 

Let £ n be the generator of ry n . Then for < s < t 

= -(ntt»)i£ (t _ a) «,(</>)) (a, x) + K(t-,) (^<^)) ( a ' *) 
= -(n^) J R" (t _ s) « s (^))(a,x) 

+it£ (t _ s) (n£ n lCs (1 - e"*) + nA((l - < s 0) - F n (l - 0(0, •)))(«, *)) 

+K (i - S ) Q nS drn£"(^ n V r ((l-<(0))-^n(l-<WO,-))) (a,x) 

= ^ (t _ s) n (-A(C(^)(0, •) - + A((l - - F n (l - 00(0, •)))(«, a:)) 

= -K (i _ s) (n^ n (0))(a,x), 

Integrating both sides with respect to s from to t, we obtain that 
(6.9) < t (0)(a,x)=K t (l-e^)(a,x)- f dsK (t _ s) (n^ n KJ))(a,x). 



o 

If 7r ni , is a Poisson random measure with intensity nis, then 

Therefore if we set ^"(0) = nii,JJ ; (^). From (j6.9|) . it is easy to see that Uf(<p) 
satisfies ([631) • D 

For any / : R + x R -» R, , we let || / ||oo= sup( 0)a! ) 6a+xR | f(a,x) \ . 
With a little abuse of notation we shall let || / ||oo= sup,j. gR | f(x) | when 
/ : R -> R as well. 



Proposition 6.3. Lei e > 0. ^ G C ; + (R + x R d ) and u^(<f>) be as in Propo- 
sition WM and ut((j>) be as in Theorem \2.4\ Then for t > e, 

(6.10) sup | u?(0)(a, x) - «t(0)(a?) H 

(a,z)eR+xR 
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Proof: For any real u E R, define, e n (u) = Xn 2 (F n (l - a) - (1 - »)) - u 2 . 
So, 



nr(</))(a,x) = J R>(l-e--)(a,x)-A / d S iO t _ s) (n^„(^))(a, x) 

Jo n 



R n nt n{l - e-n)( a ,x) - J ^ (t _ s) ( £n «(0(O-))))(a, x) 
-A / d S i^ (t _ s) (<0(O,-) 2 )(a,x) 



o 



Now 

u^(0)(a,x)-«t(^)(x) 



d S ^ (t _ a) (e„«(0(O-))))(a,x) 
+ X J^ ds (u {t _ s) «u s cf>) 2 )(a, x) - R% (t _ s) ««^(0, -) 2 )(a, x)' 

i£ t (n(l - e-^))(a,x) - £/^)(x) - jf ds^ (t _ s) (e n «(0(O, 
+A|' ( l< M (( l[s tf- U >(0 1 -) 2 )(a, 3 ;) 

+A /"ds (V s (u s </>) 2 )(x) - R n n{t _ s) (u s ^f){a,x)) 



Observe that, R n is a contraction, || u n (<j)) ||oo<|| ||oo and || u.(</>) ||oo<| 
| loo for (j> E Cy(M+ x R). Therefore, we have 

\\ uU<t>) ~ Hoc < H^Kl-e-^))- |U || e n «(<K0, •)) 

+2A || Hoc / ds || u™(4>)- u s (0) Hoc 

+A^ ds || (^_ S -K (t _ s) )(u s 



) 1 1 oo 



For (j) E C/(R+ x R rf ), note that, Ut, is a strongly continuous semi-group 
implies that u s (<p) is a uniformly continuous function. So using Proposition 
6.31 the first term and the last term go to zero. By our assumption on F, 
| e n (u™((f)(0, •)) ||oo will go to zero as n — ► oo. Now using the standard 
Gronwall argument we have the result. □ 
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Proposition 6.4. Let e > 0. The processes y n are tight in D([e, oo), M(M+ x 
R)). 

Proof By Theorem 3.7.1 and Theorem 3.6.5 (Aldous Criterion) in [7] , it 
is enough to show 

(6-11) <K+^> " <XU> ^+ 0, 

where <$> 6 C Z + (M+ x M), 5 n is a sequence of positive numbers that converge 
to and r n is any stop time of the process y n with respect to the canonical 
filtration, satisfying < e < r n < T for some T < oo. 

First we note that, as (y™, 1) is a martingale, for 7 > by Chebyschev's 
inequality and Doob's maximal inequality we have 

(6.12) P«K^> > 7) < -ci || 4> Hoc E( sup (;K\i» < -c 2 || <p |U 

7 e<t<T 7 

By the strong Markov Property applied to the process y n we obtain that 
for a, (3 > 0, we have 



L n (8 n ;a,(3) = E(e^{-a{y^ +5n ,4>) - p(y? n , $)) 
= E(exp(-(yi,ul(a<l>)+m) 
= E(exp(-(yl_ e ,u^ul(a^+m)) 



Therefore 



L n (0;a,/3) - L n (5 n ; a, (3) \< 

< ||<W n H+^)-<((a + W) lU^sup^,!}) 

t<T 

< cx || <(<(a0)+/3<£)-<((a + /3)0) |U 



where is the last inequality is by Doob's maximal inequality. Now, 



<(u£>0) + /^)-<((a + m |U<|| i?" e «>0)-a0) |U + 

+02 || ||oo / da KKW) +/30) -<((<* + Hoc +<*„(», 

JO 
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where d n {(ft) = A £ da || e n «(it^(a0) + (3(ft) + e n «((a + (3)(j>)) ||oo . 
Observe that 

II RnA u1 s n { a 4>) -oc<t>) | |oo < || Rn € {us n (acl)) - R^ Sn (a(ft)) I loo 

+ || R n ne {R n nSn {a4>) - a<f>) |U 

< || U n 5n {a4>) - RZsje 2 <P) Hoc + || K(e +5n )M) " Ke(**<t>) Hoc 

< || iC>(l - e *) - acf>) |U + / da || fl» (n 2 *(^)) |U 

+ 11 K (e+<5n) (a0)-i^ e (a«^) Hoc 

< || n(l - e#) - a<P IU +<WII |& +1)+ II K(e+5 n ){u<t>) ~ ^e(^) II 
= e n (^>) 

Consequently, 

|| + /3<A) - t£((r + !!«,< e n (<£) + d n ((f>) 

+c 2 || ||oo / da || (a$) + 13(f)) - + s)</>) ||oo . 



JO 

By Proposition 16. 1\ e n ((ft) — > and d n ((ft) — > by our assumption F n . Hence 
by a standard Gronwall argument we have that, 

(6.13) | L n (0;s,r) - L n (5 n ;s,r) |-> 

By (16+2]) . {(3^ n , 0);n = 1, 2, . . .} is tight in R + . Take an arbitrary sub- 
sequence. Then there is a further subsequence of it indexed by {n^k = 
1,2,...} such that (y2 k , (ft) converges in distribution to some random limit 
b. Thus we get 

(X£ (<f>),y% h ((ft)) (b, b) as k^ oo. 

But (|6.13p implies that 

O^M^^M) -± (6,6) as fc - oo. 

This implies that Q£* +<5n ,0) - (3 ; ™ fc fe ,<A) as k -» oo. So ([6TTD 
holds and the proof is complete. □ 

Proof of Theorem 12. 41 Proposition 16. 31 shows that the log-Laplace func- 
tionals of the process y™ converge to 3^ for every t > e. Proposition 16.41 
implies tightness for the processes. As the solution to (|2.8p is unique, we 
are done. □ 
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